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Spinon confinement in the one-dimensional Ising-like antiferromagnet SrCo2V208 
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For quasi-one dimensional quantum spin systems theory predicts the occurrence of a confinement 
of spinon excitation due to interchain couplings. Here we investigate the system SrCo2V208, a 
realization of the weakly-coupled Ising-like XXZ antiferromagnetic chains, by terahertz spectroscopy 
with and without applied magnetic field. At low temperatures a series of excitations is observed, 
which split in a Zeeman-like fashion in an applied magnetic field. These magnetic excitations are 
identified as the theoretically predicted spinon-pair excitations. Using a one dimensional Schrodinger 
equation with a linear confinement potential imposed by weak interchain couplings, the hierarchy 
of the confined spinons can be fully described. 

PACS numbers: 75.30.-m,75.10.Pq,71.70.Ej 


Confinement of quasiparticle, a phenomenon mostly 
know from particle physics, [l[ can be also realized in 
condensed matter. [H, Q Due to the strong interaction, 
quarks in a meson are asymptotically free at shorter 
distance and never exist as individual particles. The 
counterpart in condensed matter, magnetic quasipar¬ 
ticle (spinon), exhibits similar confinement behaviors. 
The spinon confinement can be illustrated in a paradig¬ 
matic one-dimensional system, Ising-like XXZ spin-1/2 
chain, where detailed theoretical descriptions have been 
achieved. In a single Ising chain, two spinons (domain 
walls) are created by a spin flip. The spinons, each car¬ 
rying spin-1/2, can move freely along the chain by sub¬ 
sequent spin flips without cost of energy. This results 
in a highly degenerate first excited state with energy J, 
the intrachain nearest-neighbor exchange interaction. A 
finite transverse component of the exchange interaction 
lifts the degeneracy of the first excited states, leading to 
an excitation continuum composed of two spinons which 
propagate independently. Q In the presence of weak in¬ 
terchain interactions, two spinons are bound because the 
separation between them will frustrate interchain inter¬ 
actions [Fig. [Il[b)-(d)]. The interchain interaction plays 
the role of an attractive potential, proportional to the 
distance, confining the spinons into pairs. This leads to 
a quantization of the excitation continuum into discrete 
levels of spinon bound states [Fig. [T][e)]. 

In this letter, we demonstrate experimentally such 
spinon confinement in the weakly coupled spin-1/2 Ising- 
like antiferromagnetic chain SrCo 2 V 208 by high reso¬ 
lution terahertz spectroscopy. SrCo 2 V 208 crystalizes 
in a tetragonal structure with the space group lAicd. 
The screw spin chains in SrCo2V208, consisting of edge¬ 
sharing CoOe octahedra, propagate along the crystalline 
c axis, as shown in the Fig. m and Fig. ^b). HQ Weak 
interchain couplings stabilize long-range collinear antifer¬ 
romagnetic order below Tat = 5 K. [HI Spin-orbit coupling 


entangles the orbital and spin degrees of freedom of Co^+, 
making the total angular momentum a conserved quan¬ 
tity. In the distorted octahedral environment, crystal 
field splitting is very strong that the lowest-lying Kramers 
doublet of Co^+ is well separated from the higher-lying 
spin-orbit quartet and sextet, thus the magnetic moment 
of Co^+ can be described by a highly anisotropic total 
angular moment with pseudospin S = 1/2.0 Significant 
Ising-like anisotropy along the c axis [Fig. [Il[d)] is also 
achieved by cooperative effects of the crystal field split¬ 
ting and spin-orbit coupling. 

Intrachain exchange interactions between the pseu¬ 
dospins in SrCo 2 V 2 08 can be modeled by the Ising-like 
XXZ Hamiltonian [il 

%,)], ( 1 ) 

where J > 0 is the nearest-neighbor antiferromagnetic 
exchange interaction, and 0 < e < I takes into account 
the Ising-like anisotropy arising from distortions of the 
CoOe octahedra. 

High-quality single crystals of SrCo 2 V 2 08 were grown 
in Helmholtz-Zentrum Berlin (HZB) using the floating- 
zone method as described in Ref. [HI. Time-domain 
terahertz (THz) transmission measurements were per¬ 
formed in the spectral range 0.1 — 3.3 THz using a TPS 
Spectra 3000 spectrometer (TeraView Ltd.). References 
(empty apertures) and samples were measured as func¬ 
tion of temperatures. Power spectra were obtained by 
doing Fourier transformation of the time-domain sig¬ 
nals. Magnetic field dependent transmission experi¬ 
ments were performed in Voigt geometry with Back¬ 
ward Wave Oscillators covering frequencies from 0.32 to 
1.4 THz and a magneto-optical cryostat (Oxford Instru- 
ments/Spectromag) with applied magnetic fields up to 
3 T. Single crystals with typical size 4 x 4 x 0.5 mm^ 
were aligned with respect to the crystalline a and c axes. 
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FIG. 1: (a) Screw chain consisting of CoOe octahedra running 
along the crystalline c axis in SrCo2V208. (b) Projection of 
left- and right-handed screw chains in the tetragonal ab plane. 
Interchain exchange interaction is dominated by the antiferro¬ 
magnetic interaction J' between the nearest-neighbor Co ions 
in the ab plane, which are from the neighboring chains with 
the same chirality, (c) Shaded area indicates the excitation 
continuum of the spinon-pair quasiparticles corresponding to 
the Ising-like XXZ Hamiltonian in Eq. O-i The spinon-pairs 
follow a quadratic dispersion relation along the chain direction 
close to the reciprocal F point {Qz = 0). ii (d) Magnetic 
excitations of spinon-pairs with total pseudospin S = 1 cor¬ 
responding to an odd number of pseudospin-1/2 flips. The 
collinear antiferromagnetic order along the c axis is stabilized 
below 5 K.0 (e) Quantized spinon-pair excitation levels due 
to linear confinement imposed by interchain exchange inter¬ 
action. 


In Fig. [2fa) we show the electric-field amplitude of 
the transmitted electromagnetic wave through a refer¬ 
ence and a SrCo 2 V 208 sample as a function of time delay 
in the magnetically ordered phase at 3.4 K. The incident 
light propagating along the crystalline a axis {k || a) is 
polarized with the ac electric field parallel to the c axis 
(e^ II c) and the ac magnetic field parallel to the a axis 
II c). The transmission spectrum in the frequency 
domain [Fig. EJb)] is obtained by Fourier transformation 
of the time-domain signal. From the power spectrum 
measured below Tyr, one can notice an absorption line at 
about 0.35 THz marked as Ei. This absorption line is not 
observed above Tat- Figure [2][c) shows the power spec¬ 
trum at 6.5 K right above divided by the one mea¬ 
sured at 3.4 K. One can readily observe the sharp line Ei 
followed by a series of absorption lines E 2 ,Es^...^Eg with 
increasing energies. These absorption lines are not ob¬ 
served when the ac magnetic field is parallel to the spin 
direction (JaE || c). This points to the magnetic nature of 
the absorption modes. These modes can be assigned to 
magnetic excitations also based on their dependence on 


FIG. 2: (a) THz time domain signals of a reference and a 
SrGo2V208 sample measured at 3.4 K with the electric field 
II c, the magnetic field || c, and wavevector k \\ a. (b) 
Transmission spectrum below and above Tat in the THz spec¬ 
tral range obtained from Fourier transformation of the time 
domain signal. The spectrum measure above Tat is shifted up¬ 
ward for clarity. A dip corresponding to Ei is marked in the 
spectrum measured below Tn- The spectral range from 1.3 
to 2.3 THz and above 3 THz is dominated by broad phonon 
bands, (c) Ratio of transmission as obtained at 6.5 K and 
3.4 K as a function of frequency below the phonon band. Nine 
peaks can be observed and are indexed from Ei to Eg. Ex¬ 
citation energy Ej versus energy level index j = 1, 2, 3,... are 
shown in the inset. 


external magnetic field as discussed in the following. 

The transmission spectrum of monochromatic THz 
light is measured at 2 K as a function of applied mag¬ 
netic field along the chain direction, which is also the di¬ 
rection of the Ising-like spins in the magnetically ordered 
phase. 0 Typical spectra are shown in Fig. [3l[a). Several 
absorption modes can be observed at the resonance fields 
in the spectrum of each frequency. The resonance mag¬ 
netic fields shift with frequencies of the electromagnetic 
waves. The absorption frequencies versus corresponding 
resonance magnetic fields are plotted in Figure [3][c). One 
can clearly recognize a series of absorption modes, which 
are split in magnetic field, with energies increasing or 
decreasing with increasing magnetic field. All the lines 
follow a similar linear field dependence that can be de¬ 
scribed by a Zeeman interaction term Egw/asHS with a 
single value for the ^-factor = 5.5, fiB being the Bohr 
magneton, and H the applied longitudinal magnetic field. 
The total pseudospin S = 1 has been used due to the fact 
that two spinons each with S = 1/2 are involved. Plus or 
minus signs correspond to split modes in magnetic field, 
and are also used to label the modes in Fig. [3][a) and 
Fig.Elc). The zero-field energies obtained by extrapolat- 
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FIG. 3: (a) Transmission spectra of various frequencies as a function of a longitudinal magnetic field at 2 K. Absorption modes 
are marked by the corresponding number of energy level in Fig.(b). (b) Normalized transmission spectrum at zero field adapted 
from Fig. [2jc). (c) Energies of absorption modes as a function of applied longitudinal magnetic field. Their field dependence 
can be described by linear Zeeman term with a single ^-factor g\\ — 5.5 (see text), (d) Hierarchy of the excitation energies at 
zero field can be modeled by the one dimensional Schrodinger equation with linear confinement Eq. CjS are the negative 
zeros of the Airy function (see text). 


ing the linear field dependence to zero field are shown in 
Fig.[3jd). The zero-field spectrum of Fig.[2jc) is shown in 
Fig. [3fb) for comparison. The agreement in the energies 
for all observed excitations is clearly documented. 


The observed magnetic excitations below T/v exhibit 
clear hierarchical features: The intensity of the lines de¬ 
creases with increasing energy and the energy difference 
between neighboring excitations also decreases with in¬ 
creasing energy [see Fig. me) and inset]. These are ex¬ 
actly the features predicted by theory for spinon bound 
states due to confinement. 0 The excitations are indeed 
the spinon-pair excitations with total pseudospin S = 1 
as confirmed by their field dependence. Every spinon- 
pair excitation is doubly degenerate in zero field, and the 
degeneracy is lifted in a finite longitudinal external mag¬ 
netic field [Fig. [3l[b)(c)]. Energies of the non-degenerate 
spinon-pair excitations increase or decrease depending on 
whether the total pseudospin is parallel or antiparallel to 
the external field, respectively. 


The S = 1 spinon-pair bound-states correspond to an 
odd number of pseudo^in flips between the spinons as il¬ 
lustrated in Fig.[T][d).y, UJj Such bound states have also 
been observed in the related compound BaCo9V90«.|12| 
The other series of spinon-pair excitations with = 0 ob¬ 
served in BaCo 2 V 2 08 by neutron scattering experiments 


[l^ are not observed in these THz measurements [see 
Fig. He)]. This confirms the magnetic dipole active na¬ 
ture of the S = 1 spinon-pair excitations. For a spinon 
pair, the relative motion of one spinon with respect to 
the other along the chain can be described by the one¬ 
dimensional Schrodinger equation 

-j^ + X\z\^={E-2Eo)^ ( 2 ) 

with the linear potential A|z| taking into account the con¬ 
finement imposed by the interchain antiferromagnetic ex¬ 
change interaction J' [Fig.Hb)], i.e., A = 2J'{Sz)‘^!c. c is 
the lattice constant along the chain direction [Fig. Ha)]. 
A quadratic dispersion relation at the F point has been 
assumed, following the calculations of the two-spinon dy¬ 
namic structure factor [Fig. Hb)].^ Since the photons 
have a linear dispersion relation with the speed of light 
involved, THz spectroscopy probes only the reciprocal 
space that is very close to the F point. 

The solution of Eq. (|2l) g ives the eigenenergies of the 
spinon-pair bound states |14| 

Ej=2Eo + Q\^/^ ( — \ j = 1,2,3,... (3) 

The sequence of the bound states is specified by the pref- 
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actors Q that are the negative zeros of the Airy function 
Ai{-Cj) = 0, Ci = 2.338,4.088,5.520,... (j = 1,2,3,...). 
2E'o is the energy threshold for creating a free spinon 
pair. The bound-state energies at zero field exhibit a lin¬ 
ear dependence on (^j as shown in Fig. [3fd). The solid 
line in Fig. [3{d) is obtained by fitting the experimental 
results to Eq. ©■ One can see that the energy hierarchy 
of the spinon-pair bound-states is nicely modeled with 
2Eq = 0.96 meV and = 0.33 meV. 

The energy threshold 2E'o is significantly smaller than 
the pure Ising-limit J = 60K:^5.17 meV. The value has 
been obtained by fitting the magnetic susceptibility using 
a pure Ising model. [HI, 111 The fact that 2Eq ^ J stems 
from the finite e term in Eq. O-lii The deviation of 
e from 1 is a measure of the anisotropy arising from dis¬ 
tortions of the CoOe octahedra, which is also reflected 
by the g\\ factor. If we assume a trigonal distortion with 
strength (5, mill then the induced anisotropy depends 
on the renormalized parameter 5/K with A < 0 being 
the spin-orbit coupling. The obtained value of g\\ = 5.5 
corresponds to 5/K = —0.80 and e = 0.73. 0 A smaller 
value of e = 0.41 in BaCo 2 V 2 08 indicates a larger distor¬ 
tion of the CoOe octahedra. 00 This is also indicated 
by a larger value of g\\ = 6.2 and a larger excitation gap 
of 1.7 meV in BaCooVoOs with a similar intrachain ex¬ 
change interaction. [3, [Hf [l2, 15, 0 

Spinon bound states in the quasi-one-dimensional an- 
tiferromagnets SrCo 2 V 2 08 and BaCo 2 V 208 can be clas¬ 
sified by their pseudospins which are either S = 1 or 
*5 = 0, analogous to the classification of mesons (quark- 
antiquark bound states) by their isospins (either 1 or 0). 
This is in clear distinction to the quasi-one-dimensional 
ferromagnet CoNb 206 where pseudospin is not a con¬ 
served quantity, since more spins are flipped at higher 
energies. Even at a certain energy level, spin blocks (do¬ 
mains) with different lengths in CoNb 206 tend to flip 
with comparable probablities. 0 

In the strong interaction scenario, the strings between 
quarks can snap and heavy hadron particles are expected 
to decay into lighter ones, when creating heavier parti¬ 
cles is energetically less favorable. Analogous to this sce¬ 
nario, spinon bound states can be expected only up to 
an energy of , above which an energy continuum 

would appear, as reported in CoNb206.[l8| This con¬ 
tinuum is surprisingly not observed in SrCo 2 V 2 08 or in 
BaCo 2 V 2 08 up to the energy much higher than Ei-\- Ei. 
The difference is possibly due to the different strengths of 
spinon-spinon interactions in the two systems. The cor¬ 
responding parameter \/y/Jl in SrCo 2 V 208 is larger by 
one order of magnitude than in CoNb206. The Ei -f Ei 
states with pseudospin S = 2 and any other bound states 
created by higher-order confinement processes with larger 
pseudospins, such as bound states of Ei states in neigh¬ 
boring chains, 0 can be excluded in SrCo 2 V 208 as doc¬ 
umented in Eig. ©c), because a stronger field dependence 
expected for these states is absent. 


Using THz transmission spectroscopy also in external 
magnetic fields, we have investigated low-energy mag¬ 
netic excitations in the quasi-one-dimensional Ising-like 
antiferromagnet SrCo2V208. Spinon-pair excitations on 
the antiferromagnetic ground state have been observed 
in this XXZ-type antiferromagnet. Spinon-pair bound 
states with entangled spin-orbit moment S = 1 are de¬ 
termined unambiguously. The hierarchy of the spinon- 
pair bound-states can be described by a one-dimensional 
Schrodinger equation with a linear confinement potential 
imposed by the interchain interaction. Energy continuum 
of the elementary spinon-pair excitations is surprisingly 
not observed in SrCo 2 V 208 . 
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